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摘 要 :基于 雷诺 平均 N-S 方程 对 脉冲 射流 作用 下 的 翼 型 非 定 常 流 场 进 行 数 值 模拟 ,采用 本 征 正 交 
分 解 (POD ) 方法 对 低频 率 、 中 等 频率 高 频率 脉冲 射流 尾 迹 中 涡 结 构 的 变化 进行 分 析 。 结 果 表 明 : 
借助 POD 方法 能 够 有 效 提取 尾 迹 中 小 尺度 涡 结 构 的 运动 状态 。POD 方法 提取 的 低 阶 模 态 主要 反 
映 了 尾 迹 中 强度 最 大 的 尾 涡 静态 分 离 结 构 , 对 应 脉冲 射流 的 主 频 成 分 ,高 阶 模 态 则 主要 反映 了 尾 滴 
内 部 或 尾 涡 之 间 的 流动 状态 ,对 应 脉冲 射流 的 高 阶 倍 频 成 分 , 尾 涡 内 部 的 相互 拉 伸 包含 频率 范围 较 
广 , 尾 涡 之 间 的 罕 引 分 解 包含 频率 较为 单一 ;低频 率 下 尾 迹 中 仅 能 提取 到 一 对 主要 尾 涡 ,高 频 座 下 
尾 迹 中 能 提取 到 多 对 强度 相当 的 尾 涡 ;低频 率 下 尾 涡 单一 且 相 互 作 用 简单 ,用 前 六 阶 模 态 即 可 表征 
流 场 结构 ,高 频率 下 尾 涡 数 量 多 且 相 互 作用 复杂 ,需要 更 高 阶 模 态 才能 表征 其 流 场 结构 。 
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Abstract:The unsteady flow field of an airfoil under the action of pulsed jet is numerically simulated 
based on the Reynolds average N-S equation. The proper orthogonal decomposition ( POD) method is used 
to analyze the variation of small-scale vortex structure in the wake of low-frequency , medium-frequency , 
and high-frequency pulsed jet. The results show that with the help of the POD method ,the motion state of 
the wake small-scale vortex structure can be effectively extracted. The low-order modes extracted by the 
POD method mainly reflect the static separation structure of the wake vortex with the greatest intensity in 
the wake , corresponding to the main frequency component of the pulsed jet. The higher-order modes mainly 
reflect the flow state inside or between the wake vortices , corresponding to the high-order frequency dou- 
bling components of the pulsed jet. The mutual stretching within the wake vortex contains a wide range of 
frequencies , and the traction decomposition between the wake vortices contains a single frequency. Only a 
pair of main wake vortices can be extracted from the wake at low frequency , and many pairs of wake vorti- 
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ces with the same intensity can be extracted from the wake at high frequency. At low frequency ,the wake 


vortex is single and the interaction is simple. The first six modes can be used to characterize the flow field 


structure. At high frequency ,the number of wake vortices is large and the interaction is complex ,so higher 


order modes are needed to characterize the flow field structure. 
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现在 飞行 器 的 发 展 历程 中 发 挥 了 重要 作用 ,但 是 复 
杂 的 舵 面 机 械 系 统 给 飞行 名 市 来 了 结构 复杂 、 维 护 
性 和 隐身 性 较 差 等 一 系列 问题 。1904 年 普 朗 特 
利用 吸 气 的 方式 控制 流体 分 离 ,主动 流动 控制 技术 
研究 从 此 开始 逐步 发 展 ” 。 

主动 流动 控制 技术 是 指 主动 控制 流体 在 物体 周 
围 或 者 内 部 流动 ,对 物体 施加 相对 应 的 升力 .阻力 或 
力 蚜 。 相 较 于 被 动 流动 控制 技术 最 明显 的 优点 是 可 
以 根据 运动 过 程 中 的 实际 条 件 , 以 最 佳 方式 进行 调 
苞 利 控制 ,实现 系统 性 能 的 改进 。 目 前 国内 外 已 
经 进行 了 大 量 关 于 主动 流动 控制 射流 的 风 洞 实验 和 
数值 仿真 。Jones 和 Engle 等 证 明了 采用 环 量 控 
制 的 机 村 可 以 获得 与 传统 高 升力 相当 ,甚至 更 高 的 
乔 汶 ,该 技术 逐 订 显现 出 缩短 起 飞 独 陆 滑 跑 距 离 、 提 
高 飞行 器 机 动 性 的 优越 性 。 兰 利 研究 中 心 ” 设计 
三 种 半 跨 模型 ,进行 了 主动 射流 控制 技术 的 跨 音 
速 实验 研究 ,人 研究 表明 后 缘 射 流 能 够 积极 改变 上 者 
面 激 波 位 置 ,减少 流动 分 离 。 北 约 科 技 组 织 的 AVT- 
239= 任 务 组 ”" 则 以 辟 身 融合 体 作为 基本 布局 ,设计 
了 物 载 主动 射流 控制 技术 的 无 人 飞行 锅 , 正 在 进行 
相关 稳定 性 和 控制 特性 人 研究。 国内 的 北京 航空 航天 
大 学 "西北 工业 大 学 “| 、 空 军工 程 大 学 “等 单 
位 也 在 定 稼 射流 参数 优化 、 增 升 减 阻 机 理 等 方面 进 
行 了 一 宇 的 实验 和 数值 模拟 研究。 

但 是 ,射流 产生 需要 从 发 动机 引 气 ,导致 额外 的 
流量 损失 。 近 年 来 有 学 者 提出 了 使 用 脉冲 射流 来 代 
蔡 定 常 射流 ,由 此 来 降低 射流 需要 的 质量 流量 和 
经 济 损失 。2006 年 在 德国 柏林 召开 的 主动 流动 控制 
ZU 总结 了 主动 流动 控制 技术 的 优势 和 存在 的 问 
题 。 会 议 指出 脉冲 射流 相 较 于 定常 射流 能 够 提供 相 
等 甚至 更 高 的 控制 能 力 ,但 是 控制 输入 却 要 小 的 多 。 

脉冲 射流 通过 周期 性 的 输入 来 调节 流 场 结 构 ， 
作用 下 的 流 场 往 往 具 有 复杂 的 流动 特征 , 尾 迹 与 外 
流 相 互 作 用 , 流 场 涡 结 构 复 杂 。 除 了 受 激励 频率 影 


响 之 外 ,加 型 绕 流 还 存在 尾 绿 涡 脱落 频率 、 边 界 层 不 
稳定 频率 、 气 流 再 附 时 涡 释 放 频 率 ”等 。 本 征 正 交 
分 解 方法 (proper orthogonal decomposition, POD ) 是 
一 种 典型 的 模 态 分 析 方 法 ,基于 主 成 分 分 析 的 原理 ， 
它 可 以 将 流 场 用 多 个 包含 流 场 主要 特征 的 本 征 模 态 
来 代 符 。 本 征 醒 态 的 分 解 按 照 能 量 大 小 的 顺序 排 
JJ ,不 同 的 模 态 能 够 反映 流 场 中 的 相关 流动 结构 , 往 
往 前 儿 阶 模仿 惑 能 反映 出 流 场 的 主要 特征 。 因 而 
POD 方法 被 广泛 应 用 于 非 定常 流 场 的 降 阶 建 模 和 多 
种 流动 问题 的 流 场 特性 分 析 “”。 流 场 特 性 分 析 已 
经 用 于 迎 角 分 离 状 态 下 的 竹 型 尾 迹 滑 分 析 、 拌 振 模 
态 分 析 ` 涡 波 流 场 分 析 OUT 等 。 

本 人 研究 基于 雷诺 平均 N-S 方程 的 数值 模拟 和 
POD 方法 ,人 研究 了 不 同 脉冲 频 座 作用 下 的 后 绿 流 场 
和 尾 迹 涡 的 变化 情况 ,对 脉冲 射流 作用 下 的 非 定 稼 
流 场 进行 分 解 , 提 取 其 主要 的 流动 模 态 ,分 析 其 作用 
规律 和 发 展 趋 势 ,对 脉冲 射流 的 理论 分 析 和 实际 应 
用 具有 一 定 的 指导 意义 。 


1 数值 计算 方法 与 验证 
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题 ,数值 模拟 方法 可 采用 二 维 雷诺 平均 Navier-Stokes 
方程 ,其 积分 形式 为 


Ef eav + f Fe ‘ndS = | Fy nds (1) 


其 中 :2 WRR 5 Fo M Fy TARAR J AR E E 
TUR TES E n D T AS TRE P IAS Eo MARK 
解 采 用 商业 软件 FLUENT ,选取 k-w SST( sheer stress 
transport ) 澳 流 模型 ,该 模型 对 于 有 较 大 逆 压 梯度 的 
边界 层 流动 分离 预测 性 能 较 好 ,采用 有 限 体 积 法 离 
散 控 制 方程 ,无 共性 通 量 采 用 三 阶 MUSCL fX, hs 
性 通 量 采 用 二 阶 迎 风格 式 离散 。 远 场 边 界 为 压力 远 
场 ,壁面 边界 采用 无 消 移 壁 面条 件 。 

本 研究 采用 文献 16 中 的 相关 模型 和 实验 数 
据 , 辟 型 相对 厚度 为 17% , 5 K c 2 240 mm ,射流 口 
mE RA-0.001c,3E ET 387 ze, K1 为 性 型 外 形 
相关 参数 。 计 算 条 件 为 :来 流 迎 角 0? ,来 流速 度 
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V, 230 m/s, 来 流 温 度 了 2293.15 上 ,雷诺 数 Re = 
500 000; 计 算 区 域 选 取 疲 型 弦 长 的 30 倍 ,网 格 训 分 
时 采用 结构 网 格 , 生 成 0 型 网 格 拓扑 结构 ,边界 层 内 
第 一 层 的 网 格 高 度 约 为 1 x10“”m, 射 流 口 附 近 第 一 
层 网 格 高 度 2 x 107^ m, PRESS — Ee RE IJ y * 均 小 
于 1 ,以 满足 黏 性 底层 的 计算 要 求 。 网 格 总 数 约 为 
58 万 ,图 2 为 六 型 网 格 划 分 和 情况。 网 格 的 收敛 性 已 
在 文献 [23 ] 得 到 验证 。 
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图 2 3879 p] fe XI A3 T DG 
Fig.2  Airfoil's computational grid 
C h FI OE eb vt Ee Pe A TE IZ cnp DC AC 
更 入 ,因此 本 研究 采用 方 波 脉冲 射流 进行 相关 研究 。 
加 其 中 脉冲 射流 在 一 个 周期 对 应 下 的 射流 速度 为 


aXiv 


Mns Le (t,t, e 
VD =i, i (2) 
jet. min te a , ) 
t 
— -open 
DC = am (3) 


式 中 :了 表示 脉冲 周期 ;tm 表示 一 个 脉冲 周期 中 射 
流 口 产 生 射 流 的 持续 时 间 ; DC 表示 占 空 比 ;V, ma K 
示 在 一 个 周期 中 的 射流 最 大 速度 ;六 .nm 表示 在 一 个 
周期 中 的 射流 最 小 速度 。 脉 冲 射 流 对 应 的 时 间 平 均 
- us [了 CDI h 
P 0 


~ dt 4 
0. 5p VS ~ 


. T 
M av. jet Ez fov, ()hdt (5) 


jet 
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前 方 来 流 静 压 为 射流 出 口 静 压 。 本 人 研究 所 用 脉冲 射 
流 的 占 空 比 DC = 7095, 射流 速度 Vau = 
113.6 m/s, 动 量 系数 C, =0.02。 本 人 研究 后 续 人 研究 均 
采用 以 上 来 流 条 件 , 仅 改变 射流 脉冲 频率 。 

目前 公开 的 轻型 高 频 射 流 实验 结果 较 少 ,因此 
选择 文献 [6 |] 中 35 Hz 脉冲 频率 作用 下 的 相关 结 
进行 验证 。 图 3 给 出 了 35 Hz 脉冲 射流 频率 下 的 时 
均 压力 系数 分 布 图 。 从 图 中 来 看 ,数值 模拟 得 到 的 
荆 型 上 下 表面 的 时 均 压 力 系 数 与 实验 结果 基本 吻 
t ,表明 本 研究 得 到 的 仿真 结果 与 实验 数据 相 接近 ， 
数值 模拟 方法 和 网 格 能 够 满足 相应 的 计算 精度 要 
求 。 图 4 给 出 了 35 Hz 脉冲 频率 下 某 一 时 刻 的 流 场 
图 ,红线 区 域 表 示 脉 冲 射 流 后 方 的 尾 迹 区 域 。 从 流 
场 图 中 看 出 , 尾 迹 区 域 的 小 尺度 涡 很 难 通过 流 线 图 
的 方式 表示 出 来 ,需要 更 精细 的 流 场 特性 分 析 方 法 。 
本 人 研究 将 借助 POD 方法 分 析 图 中 红线 所 示 尾 迹 区 
域 小 尺度 涡 结 构 的 变化 。 
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一 一 实验 数据 
—e— 仿真 数据 
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图 3 35 Hz 脉冲 射流 下 压力 系数 分 布 对 比 


Fig.3  Time-averaged pressure coefficient distribution at 35 Hz 
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图 4 尾 迹 示意 图 
Fig.4 Wake diagram 


2 POD 方法 


POD 方法 通过 将 流 场 分 解 为 奋 干 空间 正 交 模 
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态 , 并 按照 模 态 所 含 能 量 大 小 进行 排序 ,能 够 实现 对 
流 场 的 降 阶 处 理 。POD 方法 计算 的 核心 思想 就 是 寻 
RERE |V (x) j 中 的 一 组 最 优 正 交 基 1e ,ep，， 
…| ,使 得 在 1 V (x) | 中 的 元 素 在 这 组 基 上 的 投影 最 
大 , 即 求 解 以 下 最 优 解 
H e max V VULE g'e =f (6) 
P(x) 
其 中 ,(…… ) 表 示 向 量 间 的 内 积 运算 ;向 量 p 为 POD 
最 优 正 交 基 。 上 式 表 示 了 数据 集中 所 有 元 素 的 算术 
平均 。 对 于 数据 量 非常 庞大 的 非 定常 流 场 而 言 , 直 
接 求 解 上 述 最 优 问题 会 非常 困难 。 使 用 文献 [24 
提出 的 快照 技术 能 够 大 幅度 降低 计算 对 于 内 存 的 消 
耗 , 提 高 计算 效率 。 

针对 具体 的 流 场 空间 参数 进行 POD 分 析 的 过 
gU. 

设 脉冲 射流 作用 下 的 非 定常 流 场 的 空间 变量 集 
沽 "(x,y,i(i))。 其 中 x,y,i(i) 分 别 为 非 定 常 流 
场 的 横 坐 标 、 纵 坐标 和 瞬时 时 刻 。V” (x,y,t(i) ) 表 
7f 1(i) 时 刻 流 场 的 空间 变量 ,i =1,2,3,…,M。 将 
V? Cx y tC) ) 写 为 不 同时 刻 的 矩阵 格式 , 即 


n (1) (2) (n) 
Kix yt) 11a 0 [Va Vi Ut V ] 

(1) (2) (n) 

Vii Vii i Vii 

(1) (2) (n) 

Vi» Vi m Vi» 

(1) (2) (n) 

一 Vi; Vii d LN (7) 

(1) (2) (n) 

LN Vy D LEN 

E NX we q5 


利用 上 式 可 以 将 流 场 的 二 维 空间 坐标 转换 为 适 
用 于 POD 算法 的 一 维 坐标 。 首 先 将 流 场 的 空间 变 
量 分 为 时 均值 和 脉冲 值 , 即 


V(x,y,t;) - V(x ,y,t;) Vx, y, t;) (8) 


kJ ctfi V (ox, y, t) DRA IN. 个 本 征 正 交 
模 态 eo. Cc y) 及 其 对 应 时 间 系 数 w Ct; ) 乘积 的 线性 
JI, RI 


V(x,y,t;) = 人 aCi) p(y) (9) 
为 了 求解 得 到 本 征 正 交 模 态 we;,(x,y) ,借助 求 
解 特 征 值 的 方法 ,定义 相关 甜 阵 C: 
Us. 
V(x,y,t,)] 
(10) 


G - [V(x,y,t,) Vx ,y,t;) 
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求解 得 到 窃 阵 C 的 m 个 特征 值 和 及 其 对 应 的 
特征 回 量 a;(t;)。 由 特征 问 量 组 成 的 时 间 系 数 和 矩阵 
A; 为 

A; -A[a;(t;) a(t) a,(1,)] (11) 

根据 奇异 值 分解 , UV = GA4A” ,其 中 ,A = 

diag( /A,,) KS] U 2 191,02, 0,1] 经 过 归 一 化 
处 理 后 得 到 POD RZ qo; (x y) 。 

再 求 得 时 间 系 数 和 矩阵 后 也 可 以 通过 式 求 得 相应 
的 POD 其 为 


e) = gc E 0 OG) (12) 


根据 特征 值 的 大 小 可 以 对 模 态 能 量 进行 排序 ， 
则 第 7 By POD 模 态 所 含 的 流 场 能 量 占 流 场 总 能 量 比 
值 为 


E, CAL, (13) 
e 


能 量 占 比 越 高 ,包含 的 数量 集 原 信息 越 完 整 ,一 
般 来 说 ,只 需要 较 少 的 低 阶 模 态 就 能 实现 对 原始 数 
PWA ROGEN ,选用 前 大 MES BI aJ DA fiie h E 
流 场 结构 的 主要 因 系 。 


3 ”计算 结果 分 析 


3.1 MHI POD 分 解 


图 5 给 出 了 不 同 脉冲 射流 频率 下 的 升力 系数 随 
时 间 的 变化 。 低 频率 下 (5 Hz) ,升力 系数 变化 幅度 
较 大 ,此 时 脉冲 方 波 维持 峰值 射流 速度 时 间 较 长 , 升 
力 系数 基本 能 够 达到 稳定 射流 时 的 大 小 ,射流 消失 
时 ,升力 系数 呈现 断崖 式 下降 , 瞬间 丢失 邱 约 60% 
的 升力 ,随后 继续 下 降 ,直到 下 一 周期 峰值 射流 速度 
再 次 产生 。 中 等 频率 下 (50 Hz) ,升力 系数 变化 幅度 
有 所 减 小 ,射流 消失 后 ,升力 系数 的 下 降 呈 现 不 规则 
的 折线 状态 ,但 仍然 存在 断 旦 式 下 降 的 过 程 。 高 频 
Z P (200 Hz) ,升力 系数 的 变化 幅度 进一步 减 小 ,此 
时 射流 消失 后 ,升力 系数 随时 间 基 本 呈现 线性 下 降 
的 趋势 ,整体 变化 趋势 较为 规则 。 在 3/4T 时 ,在 低 
频率 和 中 等 频率 脉冲 射流 下 ,升力 系数 下 降 明显 ;而 
在 高 频率 下 ,升力 系数 下 降 并 不 明显 。 图 6 给 出 了 
不 同 脉冲 射流 频率 下 升力 系数 变化 的 FFT 分 析 图 。 
其 频率 分 布 主要 对 应 了 射流 主 频 和 其 高 阶 倍 频 ,其 
中 低频 率 对 应 的 高 阶 倍 频 分 布 更 广 ,在 流动 中 存在 
更 多 高 阶 频 率 的 成 分 。 
:应 用 力学 学 报 
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利用 POD 方法 进行 脉冲 射流 作用 下 的 辟 型 后 
绿 尾 迹 的 模 态 分 解 ,POD 分 解 的 前 30 阶 模 态 能 量 分 
布 如 图 7 所 示 , 图 7 FAS RN POD 模 态 分 解 的 阶 
数 , 纵 坐 标 表 示 各 阶 模 态 特征 值 占 总 特征 值 的 比例 ， 
即 表示 各 阶 模 态 所 占 能 量 比例 。 低 频率 下 ,一 阶 模 
态 所 占 能 量 较 高 ,前 六 阶 模 态 占据 总 能 量 的 99906 ; 


学 学 S 


(b) 中 等 频率 (50 Hz) 


图 5 不 同 射 流 频 率 下 升力 系数 的 变化 


第 39 15 
中 等 频率 下 ,一 阶 模 态 所 占 能 量 降 低 , 前 十 阶 模 态 占 
据 总 能 量 的 98. 8% ; 高 频率 下 ,一 阶 醒 态 所 占 能 量 
进一步 降低 ,前 十 阶 模 态 占据 总 能 量 的 99.4% 。 中 
等 以 及 融 频 率 下 ,一 阶 模 态 所 占 频 率 较 低 ,需要 更 多 
的 高 阶 模 态 才能 表示 脉冲 射流 的 尾 迹 变化 情况 。 


1.155 1.160 1.165 
t/s 


(c) 高 频率 (200 Hz) 


Fig.5 Lift coefficient under different Jet frequencies 
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(c) 高 频率 (200 Hz) 


图 6 不 同 射 流 频率 下 升力 系数 的 FFT 分 析 


Fig.6 FFT Analysis of lift coefficient under different jet frequencies 
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(a) 低频 率 (5 Hz) 


(b) 中 等 频率 (50 Hz) 


(c) 高 频率 (200 Hz) 


图 7 不 同 频率 下 尾 迹 的 POD 分 解 


Fig.7 POD decomposition of wakes at different jet frequencies 


3.2 尾 迹 涡 演化 及 POD 模 态 分 析 


3.2.1 IRJ Æ 
图 8 给 出 了 在 低频 率 (5 Hz) 下 ,一 个 周期 中 8 
个 瞬 态 后 绿 流 场 的 流 线 图 ,其 中 7 了 表示 脉冲 周期 。 
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在 1/8T 时 刻 , 射流 后 方 逐渐 形成 一 对 尾 涡 。 进 入 
2/8T 时 刻 , 射流 基 本 达到 稳定 状态 ,射流 后 方形 成 
完整 的 一 对 尾 涡 。 从 3/8T 到 5/8T 时 刻 , 射流 和 尾 
涡 相 对 位 置 变化 幅度 很 小 ,基本 呈现 稳定 状态 。 进 
A 6/8 T 时 刻 ,射流 消失 ,射流 后 方 的 尾 涡 被 迅速 吹 
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第 5 其 陈 款 ， 等 ， 基 于 本 征 正 交 分 解 的 脉 串 射流 尾 迹 分 析 CDimayX 


ER ,射流 对 于 竹 型 下 表面 的 阻拦 作用 消失 LARUM P € 

面 流 线 癌 上 偏 折 , 逐渐 形成 世 型 后 比 附 体 涡 。 当 进 
入 7/8T 时 刻 , 后 绿 流 场 基本 达到 无 射流 时 的 稳定 状 
fo IAEA 8/8T 时 刻 , 射流 再 次 产生 , 附 体 涡 被 压 
制 , 重 新 生成 新 的 尾 涡 。 总 体 来 看 ,在 一 个 脉冲 周期 
内 , 尾 涡 和 附 体 涡 都 经 历 了 形成 .发 展 、 消 失 .再 形成 
的 过 程 。 在 低频 率 下 ,脉冲 射流 周期 时 间 长 ,不管 是 
届 流 产生 还 是 消失 ,基本 都 能 达到 相对 应 的 稳定 
状态 。 

图 9 给 出 了 低频 率 下 前 六 阶 POD 模 态 的 速度 
云图 ,其 中 一 阶 模 态 表示 尾 迹 中 主要 的 静止 尾 涡 , 主 
要 表示 了 脉冲 射流 的 主 频 振 沪 结构 ,占据 了 整个 低 
频 脉 冲 射 流 尾 迹 的 主要 能 量 。 第 四 阶 模 态 中 出 现 尾 
PLESIR AAR, 但 直到 第 六 阶 模 态 中 尾 迹 
Z5 Th 现 相互 振荡 的 尾 涡 。 结 合 上 述 周期 中 射流 后 


Velocity(n/s): UEEIUCEEGEUUTGES 20 30 40 50 60 70 80 90100110120 
= 


er = 


图 8 5 Hz PASAM 
Fig.8 Airfoil trailing edge flow field at 5 Hz 
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方 尾 涡 的 变化 情况 ,在 低频 率 中 一 对 尾 涡 被 吹 除 并 
向 后 移动 ,由 于 下 一 对 尾 涡 再 脱落 的 时 间 周 期 较 长 ， 
尾 迹 中 仅 包含 这 一 对 尾 涡 ,这 对 尾 涡 相互 牵引 拉 伸 
并 逐渐 在 外 流 的 作用 下 向 后 移动 。 低 阶 模 态 主要 表 
示 了 该 尾 涡 的 主要 静止 状态 ,高 阶 模 态 主要 体现 了 
尾 涡 相互 拉 伸 时 的 状态 。 图 10 给 出 了 前 六 阶 POD 
模 态 时 间 系 数 的 频谱 。POD 分 解 得 到 的 模 态 往往 包 
含 多 个 模 态 ,其 中 一 阶 模 态 主 要 反映 了 脉冲 射流 的 

主 频 ,从 二 阶 模 态 开始 ,高 阶 倍 频 占 比 逐 渐 增 大 ,三 、 
四 五 阶 模 态 中 频率 构成 复杂 ,其 峰值 频率 均 为 射流 
高 阶 倍 频 ,与 升力 系数 变化 的 FFT 分 析 一 致 。 值 得 
一 提 的 是 ,在 六 阶 模 态 中 ,射流 主 频 (5 Hz) 仍 然 占据 
主要 成 分 ,这 表明 : 尾 迹 中 尾 涡 主要 以 射流 主 频 为 周 
" 隆 运 动 频率 ,其 高 阶 倍 频 主要 反映 了 尾 涡 内 部 相 


JH — 
黄 态 云图 相 一致 
牵引 拉 伸 的 状态 ,与 模 态 云 E 
EE o NU oO EN 
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图 9 5Hz PIIRRE POD EUST 
Fig.9 POD mode decomposition of wake at 5 Hz 


3.2.2. 中 等 频率 

图 11 给 出 了 在 中 等 频率 (50 Hz) 下 一 个 周期 中 
8 个 瞬 态 后 绿 流 场 的 流 线 图 。 不 同 的 是 ,由 于 射流 
持续 时 间 短 ,并 不 能 移动 到 稳定 状态 时 的 位 置 , 射 流 


后 方 的 尾 负 强度 也 较 低 ,仅仅 能 够 达到 低频 率 下 
1Z87 时 刻 的 状态 。 当 进入 6/8 T 时 刻 时 ,由 于 射流 的 
Xr ACE ,射流 口 处 的 射流 最 先 消 失 , 下 方 的 射流 由 

于 得 不 到 进一步 补充 ,逐步 被 性 型 下 表面 的 外 流 吹 
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除 ,此 时 后 缘 涡 系 复杂 , 既 有 仍 未 被 吹 除 的 射流 后 尾 


渭 , 又 有 因为 射流 口 处 射流 消失 而 先 产 生 的 附着 滴 ， 


因此 导致 中 等 频率 下 射流 消失 后 的 升力 系数 变化 复 
Zu, MEA T/8T 时 刻 时 ,此 时 后 绿 流动 状态 与 低频 


MEA 8/8T 时 刻 时 , 射流 产生 ,此 时 由 于 后 缘 流 场 
并 未 达到 无 射流 时 的 稳定 状态 , 附 体 涡 被 压制 的 过 
程 更 加 明显 。 总 体 来 看 ,中 等 频率 在 射流 消失 后 的 
后 缘 涡 系 变化 与 低频 率 差 别 较 大 。 


Vs 
率 下 6/8T 时 刻 相 似 , 王 型 下 表面 流 线 癌 上 偏 折 。 
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图 10 5 Hz F POD 模 态 时 间 功 率 谱 
Fig. 10 Spectrum of time coefficient of POD modes at 5 Hz 
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图 11 50Hz PX "m" 
Airfoil trailing edge flow field at 50 Hz 


Fig. 11 

图 12 给 出 了 中 等 频率 下 前 六 阶 POD 模 态 的 速 

度 云图 ,前 两 阶 模 态 仍 然 表 示 尾 涡 的 主要 静止 状态 ， 
从 第 三 阶 模 态 开 始 , 尾 迹 中 的 尾 涡 逐渐 呈现 相互 振 
沪 的 状态 。 这 表明 中 等 频率 下 , 尾 迹 中 存在 多 对 尾 
涡 , 但 是 只 存在 强度 最 大 的 一 对 尾 涡 ,其余 尾 涡 正 处 
于 进入 尾 迹 或 者 消散 的 阶段 。 前 两 阶 模 态 表示 了 主 
要 尾 涡 的 强度 占 总 强度 的 70 多 以 上 。 图 13 给 出 了 
前 六 阶 POD 模 态 时 间 系 数 的 频谱 。 与 低频 率 时 的 
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频谱 分 布 图 相 比 较 , 中 等 频率 下 的 频谱 分 布 网 更 加 
直观 。 一 、 二 阶 模 态 中 主 频 占 绝对 主要 成 分 ; 三、 四 
阶 模 态 中 二 倍 频 占 绝对 主要 成 分 ;五 .六 阶 模 态 中 三 
音 、 四 倍 频 占 绝 对 主要 成 分 。 
3.2.3 高 频率 

图 14 给 出 了 在 高 频 认 (200 Hz) 下 ,一 个 周期 中 
8 个 瞬 态 后 绿 流 场 的 流 线 图 。 高 频率 下 的 射流 持 终 
时 间 非 常 短暂 ,因此 导致 附着 涡 的 压制 较为 明显 ,而 
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射流 尾 涡 的 生成 则 需要 更 长 的 时 间 矿 度 , 直到 5/8T 
时 刻 ,依附 在 射流 后 方 的 尾 滴 尚未 完全 生成 。 当 进 
A 6/8 T 时 刻 ,由 于 射流 的 迟 征 效应 ,此 时 射流 后 方 


的 尾 涡 才 基 本 生成 ,由 于 缺乏 尾 涡 系 的 相互 牵引 , 因 
此 并 未 形成 中 等 频率 下 复杂 的 后 缘 涡 系 。 当 进入 
ERR, 


7/8T 时 刻 时 ,射流 完 下 型 下 表面 流 线 开 始 


NEN 
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向 上 偏 折 ,在 尚未 形成 完整 的 附着 涡 之 前 ,射流 重新 
产生 (8/87T 时 刻 )。 总 体 来 看 ,高 频率 下 ,对 于 附 体 
涡 的 压制 以 及 尾 涡 的 产生 被 放 慢 ,射流 消失 后 的 迟 
清 歼 应 变 明 显 。 相 较 于 中 等 频率 而 言 ,射流 消失 后 
的 复杂 尾 涡 系 消 失 , 迟 沸 效 应 使 升力 系数 断崖 式 下 
降 这 一 现象 消失 。 
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c Fig. 12 POD mode decomposition of wake at 50 Hz 
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图 13 


50 Hz F POD 模 态 时 间 功 率 谱 


Fig. 13 Spectrum of time coefficient of POD modes at 50 Hz 


图 15 给 出 了 高 频率 下 前 六 阶 POD 模 态 的 速度 
云图 ,与 低频 率 和 中 等 频率 不 同 的 是 ,高 频率 下 前 两 
阶 模 态 就 表示 了 多 对 尾 涡 的 相互 振荡 状态 , 且 强 度 
基本 相当 。 三 阶 和 四 阶 模 态 则 表示 了 更 多 尾 涡 的 相 
互 振荡 状态 ,这 其 中 包含 了 部 分 由 于 尾 涡 相互 诱导 
而 产生 分 解 的 洲 涡 ,五 阶 和 六 阶 模 态 则 呈现 了 更 为 
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复杂 的 尾 负 相互 运动 状态 。 图 16 给 出 了 前 六 阶 
POD 模 态 时 间 系 数 的 频谱 。POD 分 解 得 到 的 各 阶 模 
态 中 尾 涡 运 动 频率 的 分 离 变 得 更 加 明显 ,其 中 只 有 五 
阶 模 态 涉及 到 了 多 个 高 阶 信 频 成 分 。 这 表明 在 局 频 
味 流 的 尾 迹 中 ,一 对 尾 沉 中 相互 从 引 拉 伸 的 运动 状态 
减 弦 ,多 对 尾 涡 的 相互 牵引 分 解 的 运动 状态 增强 。 
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图 14 200 Hz FRR ARA 
Fig. 14 Airfoil trailing edge flow field at 200 Hz 
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图 15 200 Hz PIIRE POD 醒 态 分 解 
Fig. 15 POD mode decomposition of wake at 200 Hz 
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图 16 200 Hz F POD 模 态 时 间 功 率 谱 
Fig. 16 Spectrum of time coefficient of POD modes at 200 Hz 
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综 上 所 述 ,低频 射流 状态 下 脉冲 射流 周期 较 长 ， 
尾 迹 中 往往 只 存在 一 对 尾 滴 ;中 等 频率 下 存在 多 对 
尾 涡 ,但 仪 有 一 对 强度 较 大 ; e LS PEE S A JE 
涡 , 且 强度 基本 相当 。POD 分 解 得 到 的 低 阶 模 态 反 
映 了 脉冲 射流 尾 迹 中 主 频 对 应 的 流 场 结构 ,而 高 阶 
模仿 则 反映 了 高 阶 倍 频 对 应 的 流 场 结构 。 主 频 对 应 
的 尾 涡 主 要 反映 了 脉冲 射流 频率 的 影响 ,高 阶 倍 频 
对 应 的 尾 涡 不 仅 体 现 了 脉冲 射流 频率 的 影响 ,还 体 
现 了 尾 迹 中 尾 涡 内 部 或 者 尾 涡 之 间 的 相互 影响 。 


4 结 it 
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模拟 ,并 信 助 POD 方法 对 不 同 频率 脉冲 射流 作用 下 
的 旦 迹 进行 了 分 析 ,得 到 以 下 结论 。 

1) 高 频 脉冲 射流 作用 下 ,射流 消失 后 的 玉 游 效 
MAH T ,导致 升力 系数 断崖 式 下降 这 一 现 脓 消失 。 

2) POD 方法 提取 到 的 低 阶 醒 态 主要 反映 了 尾 
迹 中 强度 最 大 的 尾 涡 静 态 分 离 状 态 , 对 应 脉冲 射流 
主 频 的 流 场 结构 ,局 阶 模 态 则 主要 反映 了 尾 涡 内 部 
或 尾 油 之 间 的 相互 运动 状态 ,对 应 脉冲 射流 高 阶 倍 
频 的 流 场 结构 。 

3) 对 于 尾 迹 中 仅 含 一 对 尾 涡 的 低频 率 射 流 , 用 
较 纱 的 醒 态 即 可 表示 其 流动 状态 ,对 于 尾 迹 中 包含 
乡 对 尾 广 的 局 频率 射流 , 则 需要 更 多 的 醒 态 才能 够 
表示 其 流动 状态 。 避 频 率 射 流 作用 下 的 尾 迹 流 场 更 


MERET, 


4) BU UL P EPOD 方法 分 解 得 到 的 各 阶 模 
态 频 谱 岁 中 频率 分 离 情况 较 好 ,高 阶 梗 态 中 不 再 包 
含 主 频 成 分 , 尾 涡 之 间 分 解 的 运动 状态 增强 , FÉ AN 
部 相互 牵引 的 运动 状态 减弱 。 

5 ) 对 于 低频 射流 而 言 ,升力 系数 变化 幅度 较 大 ， 
HEW PF EAN [8] 5 8E EC: 0] Fé ib , AR XE Dor FH] 
T HA By KITARE ,对 于 高 频 射流 而 言 ,升力 系数 变 
化 幅度 较 小 , 尾 迹 中 存在 多 个 时 间 蜂 度 较 小 的 尾 油 。 
研究 表明 :次 型 大 迎 角 失 速 流 动 状态 的 尾 迹 中 同样 
存在 大 量 小 尺度 尾 涡 ”。 因 此 利用 高 频 射流 调节 
访 型 大 迎 角 流动 时 的 分 离 涡 具有 一 定 应 用 前 景 。 
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